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Although amphitheater-shaped valley heads can be cut by groundwater flows emerging 
from springs, recent geological evidence suggests that other processes may also produce 
similar features, thus confounding the interpretations of such valley heads on Earth and 
Mars. To better understand the origin of this topographic form we combine field ob- 
servations, laboratory experiments, analysis of a high-resolution topographic map, and 
mathematical theory to quantitatively characterize a class of physical phenomena that 
produce amphitheater-shaped heads. The resulting geometric growth equation accurately 
predicts the shape of decimeter-wide channels in laboratory experiments, 100-meter wide 
valleys in Florida and Idaho, and kilometer wide valleys on Mars. We find that when- 
ever the processes shaping a landscape favor the growth of sharply protruding features, 
channels develop amphitheater-shaped heads with an aspect ratio of n. 



1. Introduction 

When groundwater emerges from a spring with sufficient intensity to remove sediment, 
it carves a valley into the landscape (Dunne 1980). Over time, this "seepage erosion" 
causes the spring to migrate, resulting in an advancing valley head with a characteristic 
rounded form (Lamb et al. 2006). Proposed examples of such "seepage channels" in- 
clude centimeter-scale rills on beaches and levees (Higgins 1982; Schorghofcr et al. 2004), 
hundrcd-meter-scale valleys on Earth (Schumm et al. 1995; Abrams et al. 2009; Russell 
1902; Orange et al. 1994; Wentworth 1928; Laity & Malin 1985), and kilomctcr-scale 
valleys on Mars (Higgins 1982; Malin & Carr 1999; Sharp & Malin 1975). Although it 
has long been thought that the presence of an amphitheater-shaped head is diagnostic of 
seepage erosion (Higgins 1982; Russell 1902; Laity & Malin 1985), recent work suggests 
that overland flow can produce similar features (Lamb et al. 2006, 2008). To address 
this ambiguity, we seek a general characterization of processes that produce channels 
indistinguishable in shape from seepage channels. 

We first identify the the interface dynamics (Brower et al. 1983; Ben-Jacob et al. 1983; 
Kessler et al. 1985; Shraiman & Bcnsimon 1984; Marsili et al. 1996; Pelce 1988, 2004) 
appropriate for amphitheater-shaped valley heads formed by seepage erosion. We then 
show that the same dynamics apply in a more general setting. We find that whenever 
the processes shaping a landscape cause valleys to grow at a rate proportional to their 
curvature, they develop amphitheater-shaped heads with a precise shape. This result 
clarifies the controversy surrounding terrestrial and Martian valleys by showing that 



Petroff et al. 




30 m 


40 m 


50 m 


60 m 





Figure 1 . Examples of seepage valleys from the Florida network, (a) Topography obtained from 
a high-resolution map (Abrams et al. 2009) showing the rounded "amphitheater-shaped" valley 
head surrounding a spring (red arrow). Colors represent elevation above sea level, (b) A seepage 
valley head as viewed from the side wall. The red arrow shows the approximate position of the 
spring. Note people for scale. 



many of these features are quantitatively consistent with a class of dynamics which 
includes, but is not limited to, seepage erosion. 



2. The Florida network 

To provide a precise context for our analysis, we first focus on a well-characterized 
kilometer-scale network of seepage valleys on the Florida panhandle (Schumm et al. 1995; 
Abrams et al. 2009) (figure 1). This network is cut approximately 30 m into homogeneous, 
unconsolidated sand (Schumm et al. 1995; Abrams et al. 2009). Because the mean rainfall 
rate P is small compared to the hydraulic conductivity of the sand, nearly all water 
enters the channel through the subsurface (Schumm et al. 1995; Abrams et al. 2009). 
Furthermore, sand grains can be seen moving in streams near the heads, implying that 
the water drained by a spring is sufficient to remove sediment from the head. Finally, a 
myriad of amphitheater-shaped valley heads (n > 100) allows for predictions to be tested 
in many different conditions. 

We begin by finding the equilibrium shape of the water table in the Florida valley 
network. This shape describes how water is distributed between different heads. When the 
groundwater flux has a small vertical component (relative to the horizontal components) , 
the Dupuit approximation (Bear 1979) of hydrology relates the variations in the height 
h of the water table above an impermeable layer (Schumm et al. 1995; Abrams et al. 
2009) to the mean rainfall rate P and the hydraulic conductivity K through the Poisson 
equation 

yVV+P = 0. (2.1) 

To simplify our analysis, we define two rescaled quantities: the Poisson elevation <f> = 
(K/2P) 1 / 2 h and the Poisson flux q p — ||V</> 2 ||. The Poisson elevation is determined 
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entirely from the shape of the network and, consequently, can be measured from a map 
without knowledge of P or K. Physically, q p is the area that is drained by a small 
piece of the network per unit arc length. It is therefore a local version of the inverse 
drainage density (i.e., the basin area divided by total channel length). By construction, 
the groundwater flux q = Pq p . This measure of the flux differs from the geometric 
drainage area (Abrams et al. 2009) in that it is found from a solution of the Poisson 
equation, rather than approximated as the nearest contributing area. 

Fig. 2a shows the solution of equation (2.1) around the valley network (supplementary 
material). Because the variations in the elevation at which the water table emerges are 
small (<~ 10 m) relative to the scale of the network (~ 1000 m), we approximate the 
network boundary with an elevation contour extracted from a high resolution topographic 
map (Abrams et al. 2009) on which (f> is constant (see supplementary material). For a 
specified precipitation rate, this result predicts the flux q of water into each piece of the 
network. 

To test this model of water flow, we compared the solution of equation (2.1) to 
measurements at 82 points in the network. Given a reported mean rainfall rate of 
P = 5 x 10~ 8 m sec -1 (Abrams et al. 2009), we find good agreement between ob- 
servation and theory (Fig. 2b), indicating that equation (2.1) accurately describes the 
competition for groundwater. Additionally, we find that the water table elevation h is 
consistent with a ground penetrating radar survey (Abrams et al. 2009) of the area (see 
supplementary material). To understand how the distribution of groundwater through 
the network produces channels with amphitheater-shaped heads, we proceed to relate 
the flux of water into a valley head to the geometry of the head. 

3. Relation of flux and geometry 

For an arbitrary network, there is no simple relationship between the flux of water 
into part of the network and its local shape. As each tip competes with every other part 
of the network, one can only find the local flux by solving equation (2.1). However, as 
first identified by Dunne (Dunne 1980), valleys cut by seepage grow when sections of the 
valley which protrude outwards (high positive curvature) draw large fluxes while indented 
sections (negative curvature) of the network are shielded by the network. Motivated by 
this insight, we seek the relationship between the flux into a piece of a valley network 
and its planform curvature. Fig. 2c shows that this relationship is consistent with a 
hyperbolic dependence of the Poisson flux (and hence the water flux) on the curvature. 
Consequently, at tips, where the curvature is high, this relationship can be approximated 
by the asymptote. Thus, 

q p ~ Qk, (3-1) 

where fi is a proportionality constant related to the area drained by a single head. Thus 
we find a local relationship between the processes shaping a seepage valley, represented 
by the flux q p , and the local geometry, represented by the curvature k. We note that this 
relation may be further justified by a scaling argument (supplementary material), but 
here we merely employ it as an empirical observation. 

4. Geometric growth law 

In what follows we first ask how the flux- curvature relation (3.1) is manifested in the 
shape of a single valley head. To do so we first find the shape of a valley head that is 
consistent with the observed proportionality between groundwater flux and curvature. 
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Figure 2. The water table and associated groundwater flux in the Florida network, (a) The 
magnitude of the Poisson flux (color intensity on boundary) is the size of the area draining 
into a section of the network per unit length. It is found by solving equation (2.1) around 
the channels as approximated with an elevation contour. Flow lines are in black. The water 
discharge was measured at blue circles. The Poisson elevation and Poisson flux are proportional 
to the water table height and groundwater flux, respectively, (b) Comparison of the predicted 
discharge to measurements at 30 points in network taken in January of 2009 (blue points) and 
52 points in April of 2009 (red points). The black line indicates equality. This comparison is 
direct and requires no adjustable parameters, (c) We observe a hyperbolic relationship between 
the curvature of the valley walls and the predicted flux (red curve). In regions of high curvature 
(i.e. valley heads) the flux is proportional to curvature (dashed line). 



This derivation relies on three steps. First, equation (3.1) is converted, with an additional 
assumption, into a relationship between the rate at which a valley grows outward and 
its planform curvature. Next, we restrict our attention to valley heads that grow forward 
without changing shape. This condition imposes a geometric relationship between growth 
and orientation. Combining these, we find a relationship between curvature and orien- 
tation that uniquely specifies the shape of a valley growing forward due to groundwater 
flow. Finally, we find that our theoretical prediction is consistent both with valleys cut by 
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seepage and systems in which seepage is doubtful. This result leads us to conclude that 
seepage valleys belong to a class of systems defined by a specific relationship between 
growth and curvature which includes seepage erosion as a particular case. 

Following past work (Howard 1988; Abrams et al. 2009), we assume that the amount 
of sediment removed from the head is proportional to the flux of groundwater; and thus, 
by equation 3.1, Hk. From equation (3.1), the speed c at which a valley grows outward 
is therefore proportional to the planform valley curvature n. Setting H equal to the 
difference in elevation between the spring and the topography it is incising, the sediment 
flux is 

Hc = aVLn, (4.1) 

where a is a proportionality constant with units of velocity. Equation (4.1) states that 
the more sharply a valley wall is curved into the drainage basin, the faster it will grow. 
The growth of the channel head is therefore "curvature-driven" (Brower et al. 1983). 

This derivation of equation 4.1 marks a shift of focus from the mechanics that shape 
a seepage valley to the dynamics by which it evolves. Although the specific processes of 
groundwater flow and sediment transport have not been addressed explicitly, this gener- 
alization has two advantages. First, equation 4.1 is purely geometric and can be solved to 
provide a quantitative prediction for the shape of a valley head. Equally importantly, the 
generality of these dynamics suggests that the class of processes they describe may ex- 
tend beyond seepage valleys and thus provide a quantitative prediction for the evolution 
of a wider class of channelization phenomena. 



5. Shape of a valley head 

We restrict our attention to steady-state valley growth. When the channel grows for- 
ward at a speed cq without changing shape, the outward growth balances the growth 
forward. If 9 is the angle between the normal vector and the direction the channel is 
growing (Fig. 3), then c = cocos6*. Substituting this condition for translational growth 
into equation 4.1 relates the orientation of a point on the channel to the curvature at 
that point: 

cos0= (5.1) 

where n(0) denotes the dependence of curvature on orientation. Solving this equation 
(see appendix) for the shape of the curve with this property gives (Brower et al. 1983) 

y [x) = — log cos (tt- ) , (5.2) 

where w = iraQ/(coH) is the valley width and 9 = nx/w. The planform shape y(x) 
is shown in Fig. 3. A notable feature of this solution is that all geometric aspects of 
the channel head are set by the absolute scale of the valley (i.e. the valley width). In 
particular, it follows from equations (5.1) and (5.2) that all seepage channels have a 
characteristic aspect ratio 

- = TT, (5.3) 

r 

where r is the radius of curvature of the tip (Fig. 3). By contrast, a semi-circular valley 
head, in which w = 2r, has an aspect ratio of 2. 
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Figure 3. A balance between curvature-driven growth and translational growth sets the shape 
(eq. 5.2) of an amphitheater-shaped valley head (solid black curve). When a curve evolves due to 
curvature-driven growth, the normal velocity c is inversely proportional to the radius of the best 
fitting circle at that point. When a curve translates forward, there is a geometric relationship 
between the speed at which a point translates Co and the speed at which it grows in the normal 
direction c. This balance produces channels with a well defined width w and an aspect ratio of 

w/r = 7T. 

6. Comparison to observation 

To test these predictions, we first compare the shape of elevation contours extracted 
from 17 isolated, growing tips in the Florida network to equations (5.2) and (5.3). As 
these valley heads vary in size, a sensible comparison of their shapes requires rescaling 
each channel to the same size; we therefore non-dimensionalize each curve by its typical 
radius w/2. To remove any ambiguity in the position where the width is measured, w is 
treated as a parameter and is fit from the shape of each valley head. Fig. 4a compares all 
17 rescaled channels heads to equation (5.2). Although each individual valley head may 
deviate from the idealization, the average shape of all valley heads fits the model well. 

This correspondence between theory and observation is further demonstrated by com- 
paring the average curvature at a point to its orientation. We construct the average 
shape of the valley head by averaging the rescaled contours along the arc length. Rewrit- 
ing equation (5.1) in terms of the width, we obtain 



Plotting wk as a function of cos 9, we indeed observe this proportionality (Fig. 4b). 
Moreover the measured slope p = 3.07 ± 0.17 is consistent with the predicted pref- 
actor p = n. The proportionality relationship holds most clearly at high curvatures, 
where the approximation that flux scales with curvature is most accurate. Notably, were 
amphitheater-shaped valley heads semi-circular, then Fig. 4b would show the horizontal 
line wk(0) = 2. If valley heads were sections of an ellipse with an aspect ratio of it, the 



WK = 7TCOS0. 



(6.1) 
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Figure 4. The shape of valley heads in the field, experiments, on Earth, and on Mars are con- 
sistent with curvature-driven growth, (a) The shape of a channel produced by curvature-driven 
growth (red line) compared to the relative positions of points (blue dots) on the edge of valleys 
from the Florida network (17 elevation contours), (b) Comparison of the curvature at a point 
to the orientation (blue dots) of valleys from the Florida network. The red line is the linear 
relationship given in equation (6.1). The black dashed line corresponds to an ellipse with aspect 
ratio 7T. A semi-circular head would predict the horizontal line kw — 2. (c-d) The analogous 
plots for the experiments (25 elevation contours extracted at 3 minute intervals), (e-f) The 
analogous plots for three valleys near Box Canyon and Malad Gorge, (g-h) The analogous plots 
for 10 Martian ravines. 

data in Fig. 4b would follow the curve wk(6) = (4 + (tt 2 — 4) cos 2 8) 3 / 2 /n 2 . Viewing the 
semi-circle and ellipse as geometric null hypotheses, we conclude from visual inspection 
of Fig. 4b that we can confidently reject them in favor of equation (5.2). 

Seepage channels can also be grown in the laboratory by forcing water through a sand 
pile (Schorghofer et al. 2004; Howard 1988; Lobkovsky et at 2007). Because these channels 
grow on the time scale of minutes to hours, one can directly observe the development 
of the channels. Fig. 4c compares equation (5.2) to elevation contours extracted from 
a previous experiment (Lobkovsky et al. 2007) while the channel is growing. Once the 
contours are rescaled and averaged, the curvature again is proportional to cos 6 (Fig. 4d). 
The measured proportionality constant p — 3.07 ± 0.13, consistent with p = n. 

7. Generalizations 

The strong correspondence between equation (6.1) and the observed shapes of valley 
heads suggests that amphitheater-shaped heads take their form from curvature-driven 
growth. Because curvature-driven growth is a simple geometric growth model, it likely 



8 



Petroff et al. 



characterizes a class of physical processes (Brower et al. 1983). For example, when a 
low-viscosity fluid is pushed slowly into a viscous fluid in two-dimensions, the diffusing 
pressure field deforms the intruding fluid into an elongated form as described by the 
Saffman- Taylor instability (Saffman & Taylor 1958). When stabilized by surface tension, 
the shape of the resulting "viscous finger" is exactly that given in equation (5.2) (Ben- 
simon et al. 1986; Combescot et al. 1986). This morphology has also been related to the 
shape of dendrites (Mullins & Sekerka 1963; Kesslcr et al. 1986) and is a steady state 
solution to the deterministic Kardar-Parasi-Zhang equation (Kardar et al. 1986). 

This generality leads us to conjecture that when the growth of a valley head responds 
linearly to a diffusive flux, its dynamics at equilibrium reduce to curvature-driven growth. 
Geophysically relevant processes in which the growth may be dominated by a (possibly 
non-linear) diffusive flux include the conduction of heat, topographic diffusion (Culling 
1960), the shallow water equations (Chanson 1999), and elastic deformation (Landau 
& Lifshitz 1995). Thus, assuming appropriate boundary conditions exist, processes such 
as seasonal thawing, the relaxation of topography, overland flow, and frost heave may 
produce valleys indistinguishable in planform shape from seepage channels. 

To confirm the wide applicability of he geometric growth model, we proceed to com- 
pare equations (5.2) and (6.1) to enigmatic valleys on Earth and Mars. The origins of 
amphitheater-shaped heads from the Snake River in Idaho (Russell 1902; Lamb et al. 
2008) and the Martian valleys of Valles Marineris have been the subject of much de- 
bate (Higgins 1982; Malin & Carr 1999; Sharp & Malin 1975; Lamb et al. 2006). Fig. 4(e- 
h) shows that the shape of valley heads in both of these systems is consistent with equa- 
tions (5.2) and (6.1). Averaging the rescaled valleys and comparing the dimensionless 
curvature to the orientation, we find p = 2.92 ± 0.24 and p = 3.02 ± 0.21 for the Snake 
River and Martian features respectively. Both estimates are consistent with p — it. 

8. Conclusion 

That these valleys are consistent with the predictions of curvature-driven growth does 
not, however, necessarily imply that their growth was seepage-driven. We favor instead 
a more conservative conclusion: diffusive transport is ubiquitous and therefore so too is 
the log cos 9 form. 

Our results clarify the debate about the origin of amphitheater-shaped valley heads by 
placing them within a class of dynamical phenomena characterized by growth propor- 
tional to curvature. From this qualitative distinction we obtain a quantitative prediction: 
the valley head has a precisely defined shape with an aspect ratio of ir. Regardless of the 
specific mechanical processes that cause a particular valley head to grow, all valley heads 
that fall within this dynamical class will look alike. 
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1 Computation of the water table 



In order to find the distribution of groundwater flux into the network, we solved for the shape of 
the water table around the channels. From the main text, the Poisson elevation of the water table 
is a solution to the equation: 

V 2 2 + 1=O (SI) 

with absorbing and zero flux boundary conditions. Thus is independent of the hydraulic conduc- 
tivity K. 

The ground water flux at a point is related to the shape of the watertable through the equation 

?=fl|V/i 2 || (S2) 



from which 



Thus, from the definition of < 



^l|V^ 2 || (S3) 



P||V0 2 || (S4) 



Because is only a function of the network geometry, q is independent of K. This result also 
follows from conservation of mass. The total discharge from the network must be equal to the total 
rain that falls into the network, regardless of conductivity. K sets the slope of watertable at the 
boundary required to maintain this flux. 



2 Selection of the boundary 

We solve the equation around a boundary chosen to follow the position of springs and streams. 
To identify such a boundary, we first remove the mean slope (0.0025) of the topography. We 
then chose the 45 m elevation contour of the resulting topography as the boundary (Figure SI) 
obtained from a high resolution LIDAR map of the network (SI). This elevation was chosen as the 
approximate elevation of many springs. When the contour exits the area where the LIDAR map 
was available, we replace the missing section of the channel with an absorbing boundary condition. 
Because this approximation results in uncertainties in the flux near the missing boundary, we only 
analyze the water flux into a well contained section of the network (blue boundary in Figure SI). 
Finally, we include a zero-flux boundary condition in the south east in the approximate location of 
a drainage divide. We solve equation (SI) with these boundary conditions using a finite-element 
method (52). 



3 Comparison of the shape of the water table to the Poisson 
elevation 

Here we show that the solution of equation (SI) is consistent with field observations. We compare 
4> (Figure S2b) to a previously reported (SI) ground penetrating radar (GPR) survey of the channels 
(Figure S2c). 



S2 



Figure SI: To most closely approximate the shape of the network we use an elevation contour 
of the topography. Approximating the channels as nearly flat, we required that the water table 
intersect the channels at a constant height, which we chose as zero. This boundary is drawn in blue 
and black. Additionally, a drainage divide (red line) was included in the south east. Because our 
LIDAR map (SI) only shows two full valley networks, we only analyze the data from this portion 
of the boundary (blue line). The boundary is linearly interpolated between points spaced at 20 
m intervals on the blue boundary and points spaced by an average of 50 m on the red and black 
boundaries. 

As all heights are measured relative to the impermeable layer, we define h to be the reference 
elevation and shift h accordingly. It follows from the definition of <p that 



where hs is the elevation of the water table at the boundary. A least squares fit of the measured 
elevations to equation (S5) gives estimates P/K = 7 x 1CT 5 , h = 38 m, and h B = 38 m 
(figure S2d). Additionally taking P to be the observed mean rainfall rate of 5 x 1CT 8 m sec -1 , 
gives K — 6 x 10~ 4 m sec -1 . Each of these estimates is consistent with the analysis of Ref. (SI). 
Furthermore, the estimated permeability is consistent with the permeability of clean sand (S3). The 
elevation h of the impermeable layer may be overestimated due to uncertainties in the analysis of 
the GPR data. 




(S5) 



S3 




Figure S2: Comparison of the Poisson elevation to field observation, (a) The available ground pen- 
etrating radar survey was conducted on a portion of the southern valley network. The topographic 
map of the channels near the survey is 1400 m across, (b) We solved equation (SI) around the 
valley for the Poisson elevation, (c) The ground penetrating radar survey (SI) provided the ele- 
vation of the water table above sea level at 1 144 points around the network. The valley walls are 
represented by the elevation contours for 30 m to 45 m at 5 m intervals, (d) The measured height 
is consistent with theory. The red line indicates perfect agreement. 
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4 Comparison of contour curvature to the groundwater flux 



Because the curvature is a function of the second derivative of a curve, its estimation requires 
an accurate characterization of the channel shape. To have the highest possible accuracy in the 
estimation of curvature, we restrict the comparison between flux and curvature to a small piece of 
the network where the boundary is linearly interpolated between points separated by 5 m. 

The curvature at a point on the boundary is computed by fitting a circle to the point and its 
neighbors on both sides (Figure S3b). Given the best fitting circle, the magnitude of the curvature 
is the inverse of the radius. The curvature is negative when the center of the circle is outside the 
valley and positive when the center is inside the valley. 

To compare the curvature and flux at a point, we calculate the Poisson flux q p into each section 
of this piece of the network by solving equation (SI) between the channels (Figure S3c). We closed 
the boundary on the eastern side of the domain by attaching the extremities to the valley network 
to the east using zero-flux boundaries. To identify the characteristic dependence of the flux at a 
point on the curvature, we averaged the flux and curvature at points on the boundary with similar 
curvatures. Each point in Figure S3d represents the average flux and curvature of 50 points on the 
boundary. 



The Poisson flux is the area that drains into small segment of the network divided by the length of 
the segment. It can therefore be considered as a "local" inverse drainage density. Because all of 
the area drains into some piece of the channel, the integral of the Poisson flux is the total area of 
the basin. It follows that its mean value is the inverse drainage density. 

In what follows we ask how the Poisson flux depends on the distance d a piece of the network 
is from its drainage divide. We note that if a d has a characteristic value in a network, then we find 
a scaling of geometric flux with curvature that is consistent with observation (Fig. lc). 

A section of the network receives a large flux when it drains a large area a or when all of the 
water is forced through a small length of channel wall i. When water from a large basin (d 3> k^ 1 ) 
drains toward a point , then a ~ d 2 (Fig. S4a). Note that "~" is the symbol for "is the order of 
magnitude of" or "scales as." This area is drained into a section of channel, the length I of which 
is proportional to the planform radius of curvature, k -1 ; thus in regions of high curvature 



where f2 = md 2 is a constant of the network related to the characteristic groundwater discharge 
of a head and m ~ 1 is a proportionality constant related to the characteristic shape of a valley 
head. The flux into a point is therefore proportional to the product of variables characterizing the 
network, f2, and the local geometry of the channel, k. Equating d with the inverse drainage density 
of the network, we find d = 147 m from the analysis of the topographic map. Fitting a hyperbola 
to the data in Fig. lc, given this value of d, gives m — 1.5 ± 0.2, consistent with m ~ 1. 

In concave regions of the channel the area drained is the sum of the area outside the concavity 
and the area inside the concavity (Fig. S4b). This area a can be expressed as 



5 The Poisson flux-curvature relation 
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a = m\dK 1 + m 2 K 



-2 
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Figure S3: Identification of the relationship between the curvature of the valley walls and the 
local flux of groundwater, (a) The curvature and flux are measured between two valleys along the 
black contour, (b) The curvature at each point on the boundary is measured by fitting a circle to 
boundary, (c) The flux into each section of the network is found from the solution of equation (SI), 
(d) Comparison of the flux into each section of the network to the curvature. Geometric reasoning 
gives the asymptotic behavior (black dashed lines) of this relation when the magnitude of the 
curvature is large. 
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a 



a = 7V(f/2 



a\ — 2k 1 d $ 



1 = 7t/k 




£ = tt/k 



Figure S4: The Poisson flux is the local drainage density, (a) When a basin drains into a convex 
region (red line) the drainage density increases with curvature k. (b) When a basin drains into a 
concave interval, the drainage density decreases with curvature. 



where mi and m 2 are dimensionless numbers related to the shape of the drainage basin outside 
and inside the concavity, respectively. For example, if the concavity is a semi-circular depression 
and it drains a rectangular region, then mi = 2 and m 2 = n/2. This area is drained by a segment 
of length ~ k~ 1 giving a mean Poisson flux q p that scales as 

q p = (mid + m 2 K,~ 1 ) /m 3 , (S8) 

where m 3 is a dimensionless number related to the shape of the concavity. Fitting the data to a 
hyperbola, and again taking d = 147 m, we find mi /m 3 = 1. 52±0. 22 and m 2 /m 3 = 10.80±2.97. 
This scaling relation, in combination with the behavior at large positive k, gives the the asymptotic 
behavior of the flux-curvature relation. 



6 Derivation of the shape of the valley head 

Here we derive equation (5) of the main text. 

The balance between translation and curvature-driven growth relates the orientation to the cur- 
vature through the equation 

7TCOS# = WK. (S9) 

We first re- write the orientation of a segment in terms of the local normal n(x) to the curve and the 
direction the head is translating y. It follows from the definition of 9 that 

7in(x) ■ y = wk(x), (S10) 

Next, by describing the shape of a valley head by a curve y(x), equation (S10) becomes 

-7T d xx y 

w /-, , /o w , fr (SI 1) 



[d x yY (1 + (9^)2)3/2 • 
S7 



With the substitution g = d x y, this equation is re-expressed as an integrable, first order equation as 

wd x g + tt(1 + g 2 ) = 0. (S12) 

Integrating once, 

a = d~v = — tan 

V w 



/ 7VX s 

g = d x y = -taa[ — ). (S13) 



Integrating a second time for y gives 

y= -log cos — , (Sl4) 

7T V W / 

equivalent to equation (5) of the main text. 

Although not necessary here, it is occasionally useful to express the shape of the channel as a 
vector v parameterized by arc length s, 

, x w ( 2arctan(tanh(7rs/2u>)) \ 

v(s) = — . ( U ( I w \ ■ ( Sl5 ) 
7r \ log(sech(7rs/w)) y 

The derivative v is the unit tangent vector. 



7 Selection of valley heads 

The derivation of equation (Sl4) requires that the channel grow forward without changing shape. 
Consequently, when identifying seepage valley heads suitable for analysis, we restricted our anal- 
ysis to isolated channels. 




Figure S5: 17 isolated valley heads were chosen from the Florida network 
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7.1 Florida Network 



We select valley heads from the Florida network that are reasonably isolated and not bifurcating. 
Given such a valley, we extract an elevation contour approximately one half the distance between 
the spring and the upland flat plain. We find that the deviation in the shape of any given channel 
from equation (S14) is insensitive to the choice of elevation contour. 

Table SI: Valley heads from the Florida network. Coordinates are given with respect to UTM zone 
16R 

channel Easting (m) Northing (m) elevation (m) width (m) 



1 


696551.40 


3373949.52 


56.91 


100.66 


2 


696423.49 


3373123.32 


55.74 


111.03 


3 


694537.55 


3373068.53 


42.93 


50.47 


4 


693995.09 


3373701.11 


49.49 


49.60 


5 


694391.80 


3373813.01 


43.49 


38.72 


6 


696841.09 


3373900.80 


44.23 


28.96 


7 


698339.72 


3374200.55 


59.95 


83.53 


8 


698040.54 


3374282.69 


50.91 


48.28 


9 


697285.68 


3375011.47 


59.12 


80.20 


10 


695968.97 


3375029.24 


49.15 


56.05 


11 


696114.42 


3375019.47 


46.62 


42.59 


12 


696336.97 


3375135.23 


49.88 


56.34 


13 


696453.90 


3375233.09 


51.13 


50.54 


14 


696976.13 


3375317.38 


51.10 


46.08 


15 


694818.57 


3375532.39 


54.82 


52.04 


16 


698537.91 


3374777.58 


54.21 


70.77 


17 


697463.52 


3375108.63 


53.97 


55.17 
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10 cm 




46 minutes 82 minutes 119 minutes 



Figure S6: An elevation contour (blue lines) was extracted from the experiment every three min- 
utes from a digital elevation map (55). These three representative elevation contours from the 
beginning, middle, and end of the experiment demonstrate that the shape changed little during 
growth. 

7.2 Experiments 

The experimental apparatus used to grow seepage channels has been previously described (S4). 
The channel used in the comparison to equation (S14) grew from an initially rectangular indenta- 
tion 3 cm deep in a bed of 0.5 mm glass beads sloped at an angle of 7.8° with a pressure head of 
19.6 cm. To extract the shape of the channel, we first removed the slope of the bed by subtracting 
the elevation of each point at the beginning of the experiment. We then follow the growth of an 
elevation contour a constant depth below the surface. Because the shape of channel at the begin- 
ning of the experiment is heavily influenced by the shape of the initial indentation, we restrict our 
analysis to the shape of the contour after 45 minutes of growth. The channel grew for a total of 
119 minutes and was measured at 3 minute intervals. 
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7.3 Snake River valley heads 



To compare the form of amphitheater-shaped valley heads growing off of the Snake River in Idaho, 
we extract the valley shape from images taken from Google Earth. We select three prominent heads 
(Table S2, Figure S7); Box Canyon (S6) and two near Malad Gorge. We extract the shape of each 
of these heads by selecting points at the upper edge of the valley head. The mean spacing between 
points is 13 m. We stop selecting points when the valley turns away from the head. 



100 m 




head 1 head 2 head 3 



Figure S7: The shape of amphitheater-shaped valley heads growing off of the Snake River in Idaho 
were extracted from aerial photos of the channels. Heads 1 and 2 are near Malad Gorge. Head 3 is 
Box Canyon. 



Table S2: Valley heads near the Snake River 
channel latitude longitude width (m) 

1 42.8675° 115.6432° 190 == 

2 42.8544° 115.7045° 166 

3 42.7084° 114.9683° 132 
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7.4 Martian valley heads 



The shapes of the Martian ravines which we compared to equation (S14) were extracted from im- 
ages generated by the Themis camera on the Mars Odyssey orbiter. Channels are selected based on 
the condition that the amphitheater head was largely isolated from neighboring structures. Because 
the ravines are deeply incised into the topography, there is typically a sharp contrast between the 
ravines and the surrounding topography. We extract the shape of the ravine by selecting points 
spaced of order 100 m apart along the edge of the ravine (Table S3, Figure S8). We stop selecting 
points when the ravine intersects with a neighboring structure or when the direction of the valley 
curves away from the head. 

1 km 

head 1 head 2 head 3 head 4 

^ ^ n C] 

head 5 head 6 head 7 head 8 



head 9 head 10 



Figure S8: 10 valley heads near the Nirgal Valley, Mars. The shape of each head was extracted by 
selecting points at the edge of the valley head from images generated by the Mars Odyssey orbiter. 



Table S3: Martian valley heads 
head Themis Image latitude longitude width (m) 



1 


V06395001 


-8 


.7270° 


278.1572° 


4730 


2 


V06395001 


-8 


.7235° 


278.1557° 


2650 


3 


V09004001 


-9 


.4310° 


274.6110° 


1940 


4 


VI 1138002 


-7 


.9183° 


275.4740° 


3690 


5 


VI 1138002 


-7 


.9160° 


275.4736° 


3740 


6 


V14133002 


-9 


.5763° 


278.4435° 


2940 


7 


V14857001 


-7 


.5656° 


273.6060° 


3110 


8 


VI 6654002 


-8 


.7792° 


275.5868° 


3970 


9 


VI 6654002 


-8 


.7781° 


275.5894° 


3310 


10 


V26750003 


-8 


.0633° 


274.8977° 


3370 
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8 Stream discharge data 



8.1 Comparison of field measurements to the predicted flux 

Fig. 1 of the main text compares the solution of equation (SI) to field measurements. The instan- 
taneous discharge of a stream is measured from the the cross-sectional area a in a locally straight 
section of the channel and the surface velocity v, from which the discharge Q = av. We measure 
the surface velocity of the stream from the travel time of a small passive tracer between points 
at a fixed distance. This method may underestimate the discharge in very small streams where a 
substantial fraction of the flow may be moving through the muddy banks of the stream. 

To compare the measured discharge to the Poisson equation, we integrate the flux, q = P||V0 2 ||, 
along the section of the network upstream from the measurement assuming the reported annual 
rainfall, P = 5xlCT 8 msec~ 1 . When discharge is measured near a spring, the flux is integrated 
around the valley head. 
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January 2009 

Easting (m) Northing (m) discharge (cm 3 sec -1 ) predicted discharge (cm 3 sec -1 ) 



696905.45 


3374708.15 


11700 


4802.54 


695425.47 


3374595.35 


310 


422.30 


695333.80 


3374486.16 


1900 


1459.50 


695410.03 


3374422.51 


2100 


1403.02 


695589.53 


3374413.20 


2000 


1272.90 


695608.45 


3374439.65 


1700 


1145.17 


695602.26 


3374467.29 


4700 


2103.67 


695532.14 


3374764.77 


310 


490.24 


694045.68 


3373713.71 


710 


1708.72 


694102.24 


3373742.47 


850 


1381.31 


694110.98 


3373726.59 


850 


3520.83 


694393.38 


3373788.44 


810 


1761.83 


694515.20 


3373714.30 


2300 


2051.40 


694700.99 


3373494.69 


2900 


1831.06 


697174.63 


3373662.18 


700 


3004.36 


697622.18 


3374045.11 


10800 


4700.85 


697523.57 


3374034.52 


440 


2225.53 


696432.08 


3373937.74 


3500 


2619.48 


696353.61 


3374006.59 


3500 


2688.01 


696415.16 


3373979.53 


3600 


2545.77 


696363.79 


3373884.98 


570 


673.20 


696314.56 


3373838.46 


3100 


1132.34 


695400.74 


3373894.43 


2800 


2117.42 


695417.69 


3373884.87 


3100 


2522.08 


694429.25 


3374329.77 


700 


1145.16 


694541.01 


3374318.70 


1250 


1626.34 


694295.68 


3374320.27 


700 


950.41 


694081.94 


3374205.31 


1950 


2969.80 


693696.69 


3373094.27 


300000 


284251.54 


693575.95 


3374496.41 


100000 


148834.48 
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April 2009 

Easting (m) Northing (m) discharge (cm 3 sec -1 ) predicted discharge (cm 3 sec -1 ) 



696577.00 


3375064.00 


26688 


17558.21 


696515.00 


3375060.00 


31324 


25683.40 


696526.00 


3375081.00 


10878 


7960.94 


696378.00 


3375075.00 


3155 


2196.50 


696374.00 


3375047.00 


31168 


28808.33 


696312.00 


3374949.00 


44714 


33750.34 


693684.21 


3374490.78 


181142 


148233.29 


693857.70 


3374480.17 


134261 


134367.19 


694237.23 


3374550.07 


159510 


123111.22 


694371.93 


3374575.91 


96597 


120688.53 


694445.52 


3374574.44 


123230 


115351.89 


694706.00 


3374606.66 


142841 


111689.91 


694808.19 


3374666.05 


133061 


103032.24 


694815.26 


3374674.12 


24251 


14558.59 


695449.36 


3374792.20 


70782 


70543.57 


695317.14 


3374776.63 


115771 


82714.80 


695400.81 


3374783.02 


18354 


10476.11 


695613.16 


3374808.59 


46630 


69195.87 


695756.20 


3374863.59 


11422 


10024.07 


695787.02 


3374851.22 


81630 


57339.57 


695914.95 


3374827.10 


24757 


41590.16 


695922.74 


3374822.92 


31480 


15071.30 


696011.72 


3374871.04 


6903 


2472.52 


696019.12 


3374873.43 


52090 


38588.88 


696127.23 


3374876.96 


44644 


36223.61 


696267.06 


3374905.73 


51745 


34747.90 


696335.93 


3374970.34 


52171 


29296.08 


696577.00 


3375064.00 


26688 


17558.21 


696515.00 


3375060.00 


31324 


25683.40 


696526.00 


3375081.00 


10878 


7960.94 


696378.00 


3375075.00 


3155 


2196.50 


696374.00 


3375047.00 


31168 


28808.33 


696346.42 


3374960.48 


8141 


4287.83 


696916.41 


3374703.37 


2704 


2773.16 


696913.37 


3374697.05 


2131 


1230.01 


695406.38 


3373894.53 


6791 


2117.42 


695284.53 


3373820.41 


4975 


6209.90 


695268.73 


3373828.01 


12171 


7245.29 


695207.06 


3373539.81 


28499 


16435.07 


695163.07 


3373472.73 


285299 


214562.12 


695825.64 


3373844.08 


20009 


6515.55 



Easting (m) Northing (m) discharge (cm 3 sec x ) predicted discharge (cm 3 sec x ) 



695818.15 
695829.40 
695870.78 
695873.99 
694804.18 
694811.43 
694864.00 
694853.04 
694999.62 
695043.10 
695043.00 
695410.00 
697528.00 
695529.00 
695437.00 
695434.00 



3373874.10 
3373872.17 
3373925.80 
3373937.31 
3374918.55 
3374929.34 
3374985.35 
3375015.85 
3375057.93 
3375092.62 
3375070.20 
3373885.00 
3374024.00 
3374749.00 
3374602.00 
3374600.00 



3098 
6847 
1292 
7298 
4777 
15554 
9906 
6866 
11789 
3376 
5248 
4173 
995 
685 
263 
10759 



1747.84 

4724.56 

1436.14 

2562.84 

2926.35 

10799.79 

7942.65 

2769.17 

6421.11 

2175.13 

3776.13 

2590.45 

2400.35 

490.24 

422.30 

9777.40 
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